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ABSTRACT 

Recent studies of active galactic nuclei (AGN) found a statistical inverse linear scaling be¬ 
tween the X-ray normalized excess variance op ms (variability amplitude) and the black hole mass 
spanning over Mbh = 10 6 — 10 9 M 0 . Being suggested to have a small scatter, this scaling rela¬ 
tion may provide a novel method to estimate the black hole mass of AGN. However, a question 
arises as to whether this relation can be extended to the low-mass regime below ~ 10 6 Mq. If 
confirmed, it would provide an efficient tool to search for AGN with low-mass black holes using 
X-ray variability. This paper presents a study of the X-ray excess variances for a sample of AGN 
with black hole masses in the range of 10 5 — 10 6 M 0 observed with XMM-Newton and ROSAT , 
including data both from the archives and from newly preformed observations. It is found that 
the relation is no longer a simple extrapolation of the linear scaling; instead, the relation starts 
to fl atten at ~ 10 6 M 0 toward lower masses. Our result is consistent with the recent finding of 
Ludlam et al. ( 20151) . Such a flattening of the Mbh — of ms relation is actually expected from the 
shape of the power spectrum density of AGN, whose break frequency is inversely scaled with the 
mass of black holes. 


Subject headings: galaxies: active — galaxies: nuclei — X-rays: galaxies 


1. INTRODUCTION 


Rapid X-ray variability is one of the ba¬ 
sic observational c haracteristics of active galac¬ 
tic n u clei (AGN) ( McH ardvl Il985 : Grandi et al 


1992 : iMushotzkv et al. 19931) . It is a useful 
tool to study black holes (BH) and the cen¬ 
tral engine of AGN, since the X-ray emission 
is thought to originate from the innermost re¬ 
gion of an accretion flow around the BH. One 
commonly used method to characterize the vari¬ 
ability is the power spectrum density (PSD) 
analysis, which quantifies the amount of vari¬ 
ability power as a function of temporal fre- 
quency ( Green et~aDl993l : Lawrence fe Papadakis 
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19931) . The PSD of AGN has been found to 
be_well descr i bed by a bro ken p ower- law (e.g. 
Papadakis fe McHardvl 1995 : Edelson fe Nandra 


19991 lUttlev et al.l 21)0 2t ~ Markowitz et al.1 20031 : 

Vaughan et al.1 2003bl) . The break frequency is 
found to be inversely scaled with the black hole 
mass in a linear way, with a possible dependence 
on the scaled accretion ra te m (in units of the 
Eddington accret ion rate ) (iMcHardv et al.1 20061 


Go nzale z-Martin & Vaughan 2012|). These results 
are remarkable in the sense that AGN show simi¬ 
lar X-ray variability properties to black hole X-ray 
binaries (BHXB), indicating th at AGN are scale d- 
up versions of BHXB (see also IZhou et al. 2015 ). 


However, reliable PSD analyses require well 
sampled, high quality X-ray data of time se¬ 
ries, which are generally hard to obtain for large 
samples of AGN for the current X-ray obser¬ 
vatories. Instead, an easier-to-calculate quan¬ 
tity, the “normalized excess variance” (e.g. 
iNandra et al.1 [l997 ; [Turner et al.l 1 19991) , is com¬ 
monly used to quantify the X-ray variability am- 
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plitude. Early studies revealed correlations be¬ 
tween the excess variance and various param¬ 
eters of AGN, such as X-ray luminosity, spec¬ 
tral in dex, and the FWHM of the H/3 line 
' e.g. Nandra et al. 1997 ; George et al.1 200d: 


Turner et al. 1999t Markowitz fe; Edelson l200lh ' 

Later work suggested that these correlations are 
in fact by-products of a more fundamental relation 
with black hole mass: the M 


(iLu fc Yul 2 OOIL Bian fe ZhaoT 20031 : Papadakis 


h — uL relation 


2004)). In fact, this relation conforms to the scal- 


BH mass ( 

Papadakis 

O 

O 

CM 

O’Neill et al. 

20051 

Ponti et al. 

2012|). O’Neill et al. (2005li confirmed 


the anti-correlation between the excess variance 
and Mbh wi th a lar ge AGN sample observed with 
ASCA. Zhou et ahl ( 2010h obtained a tight cor¬ 
relation, using high quality XMM-Newton light 
curves of AGN whose black hole masses were 
measured with the reverberation mapping tech¬ 
nique. The intrinsic dispersion of the relationship 
(~ 0.2 dex) is comparable to that of the relation 
between Mbh and stellar velocity d ispersion for 
galactic bulge ( Tremaine et al. 2002ll . By making 
use of a large sample of 161 AGN observed with 
XMM-Newt o n for at least 10 ks for each object, 
Ponti et al.l ( 2012h (henceforth P12) reaffirmed 
this relationship (using the excess variance calcu¬ 
lated on various timescales of 10 ks, 20 ks, 40 ks 
and 80ks), but found only a weak dependence 
on the accretion rate. The significant correlation 
with small scatters (0.4 dex for the reverberation 
mapping samp l e, 0. 7 dex for the CAIXAvar sam¬ 
ple, see Ponti et al ■ 203 for details) suggested 
that it may provide similarly or even more accu¬ 
rate black hole mass estimation compared to the 
method based on t he single epoch optical spectr a 
( Kaspi et al. 2000t Vestergaard &; Peterson 2006h . 


Moreover, unlike the commonly used virial method 
that is susceptib l e to t he orientation effect of AGN 
( Gollin fc Kawaguchi! 2004 1. X-ray variability can 
be considered as inclination-independent. 

However, the previous studies were based on 
AGN samples with mostly supermassive black hole 
of Mbh > 10 6 M 0 , and little is known about the 
relation for AGN with Mbh < 10 6 M 0 . The of ms 
of a few AGN with Mbh ~ 10 6 M 0 were presented 
and found to have the larges t cr? ma values among 
AGN over a la rg e Mrh range (iMiniutti et al.l2009t 
Ai et al. 201ll f. IPonti et all ( 201 2li proposed that 


the Mbh — relation may show a deviation 
from the linear relation in the low-mass regime; 
however, the data is too sparse to draw a firm con¬ 
clusion. Thus, the question remains unanswered 
as to whether this relation can be extended to 

M bh < 10 6 Mq 0. 

The answer to this question is important in at 
least two aspects. First, if the answer is “yes”, 
the relation would provide a valuable method to 
find the so-called low-mass AGN with Mbh 
10 6 M 0 , or sometimes referred to as AGN with 
intermediate-mass black holes (IMBH). This is of 
particular interest since in these AGN the com¬ 
monly used virial method involving optical broad 
emission lines becomes difficult in practice due 
to the faint AGN (broad line) luminosities that 
are outshined by the host galaxy starlight. In 
fact, there has been a few attempts in practice 
by adopting such an assumption. By extrapolat¬ 
ing the Mbh — o')r elation to below 10 6 M 0 , 
iKanhzasa et al. ( 2012 1 selected a sample of 15 
candidate low-mass AGN using the X-ray excess 
variance. Second, the Mbh — of ms relation must 
have a cutoff somewhere, otherwise the variabil¬ 
ity would have become unrealistically large for 
ultra-luminous X-ray sources (ULX) and BHXB, 
which js not seen in o bserv a tions, however (e.g. 
Gonzalez-Martfn et al. 2011 : [Zhoul 12015 ). From 
a theoretical perspective, the Mbh — of ms rela¬ 
tion is actually a manifestation of the inverse scal¬ 
ing of the break frequency of the AGN PSD with 
the black hole mass, since the excess variance is 
the integral of the PSD over frequency domain 
( van der Klis 19891 1997 1. In fact, a break of the 
Mbh — of ms relation at low BH masses had been 
predicted based o n the current understanding of 


the PSD of AG N (IPapadakisl 2004 : O’Neill et al 
120051 : IPonti et ah 2012T) . and the exact break mass 


mass 

(the mass at which the break occurs) depends on 
the shape of the PSD. Therefore, a study of the 
Mbh — ^rms relation in the low-mass regime may 
provide a constraint on the break mass, as well as 
the shape of AGN PSD. 

In this paper, we study the Mbh — rela¬ 
tion in the Mbh = 10 5 — 10 6 M 0 range, using an 
optically selected low-mass AGN sample from our 


1 While we were w riting this paper, a new paper 
iLudlam et al.|[2015l l appeared very recently, which carried 
out a similar study and achieved a similar result as ours in 
this work. 
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previous work (Dong et alj 12012). We use both 
new observations and archival data obtained with 
XMM-Newton, as well as archival data of ROSAT. 
This paper is organized as follows. The introduc¬ 
tion of our sample and data reduction are pre¬ 
sented in section 2. In section 3, the excess vari¬ 
ance is introduced, as well as the PSD models of 
AGN concerned in this study. The results and dis¬ 
cussion are presented in section 4. Throughout the 
paper, a cosmology with Hq = 70kms _1 Mpc^ 1 , 
D m = 0.3 and Da = 0.7 is adopted. 


2. SAMPLE, OBSERVATION AND DATA 
REDUCTION 


There are over 300 low-mass (type 1) AGN with 
Mbh < 10 6 known so far. The largest is 
from our work ( Dong et al. 20121) . which was se¬ 
lected homogeneously from the Sloan Digital Sky 
Survey (DR4), and comprises 309 objects with 
Mbh < 2 x 10 6 Mg. The BH masses were es¬ 
timated from the luminosity and the width of 
the br oad Ha line, u sing the vir ial mass formal¬ 
ism of Greene fe Hoi ( 20051120071) . One feature of 
this sample is the accurate measurements of the 
AGN spectral parameters, and hence the black 
hole masses and Eddington ratios. Compared to 
previous samples, this sample is more complete as 
it includes more objects with low Eddi ngto n ratios 
down to Lboi/AEdd ~ 0.01 (see also lYuan et al 


20141) . We compile a working sample of low-mass 


AGN with usable X-ray data from this parent sam¬ 
ple. 


We search for X-ray observations from both the 
X'MM-Newton and ROSAT PSPC data archives 
to maximize the sample size. For X'MM-Newton 
observations the 3XMM-DR4 catalogue is used. 
We also add new observations of three objects 
from our programme to study low-mass AGN 
with XMM-Newton (proposal ID: 074422, PI: W. 
Yuan). We consider only observations with ex¬ 
posure time longer than 10 ks. There are 26 
XMM-Newton observations for 16 objects and 6 
ROSAT observations for 6 objects found in the 
archives. The new observations of J0914+0853, 
J1347+4743, and ,11153+4612 were performed by 
XMM-Newton at faint imaging mode on November 
1st, 22nd, and December 4th 2014 with an expo¬ 
sure time of 36 ks, 31 ks, and 14 ks, respectively. 

The X-ray data are retrieved from the X'MM- 


Newton and ROSAT data archives. We follow the 
standard procedure for data reduction and anal¬ 
ysis. For the XMM-Newton observations, we use 
the data from the EPIC PN camera only, which 
have the highest signal-to-noise. Light curves are 
extracted from observation data files (ODFs) by 
using the XMM-Newton Science Analysis System 
(SAS) version 12.0.1. Events in the periods of 
high flaring backgrounds are filtered out @. Obser¬ 
vations with cleaned exposure time shorter than 
10 ks are also excluded. Typical source extraction 
regions are circles with a 40 arcsec radius. Only 
good events (single and double pixel events, i.e. 
PATTERN < 4) are used for the PN data. Back¬ 
ground light curves are extracted from source-free 
circles with the same radius. Finally the SAS task 
EPICLCCORR is applied to make corrections for 
each of the XMM-Newton light curves. The energy 
band 0.2-10 keV is used. The time bins of the light 
curves are chosen to be 250 s, which are the same 
as in P12 for easy comparison. As demonstration, 
Figure |T| shows some of the typical light curves 
(panels 1-4). 

For ROSAT PSPC observations, the XSELECT 
package is used to extract source counts and light 
curves. Typical source extraction regions are cir¬ 
cles of 50 arcsec radius, and the same aperture is 
used to extract background light curves. The en¬ 
ergy band for ROSAT observations is 0.1-2.4keV. 
The time bins are also set to be 250 s for the same 
reason. Examples of the ROSAT light curves are 
also shown in Figure [T| (panel 5). 

All the light curves exhibit significant variabil¬ 
ity on short timescales, as shown in the figure. For 
a few sources, due to short observational inter¬ 
val and relatively low signal-to-noise ratio (S/N), 
the intrinsic variability is overwhelmed by ran¬ 
dom fluctuations because of the large statistical 
uncertainties (the statistical uncertainty is larger 
than the source variability; see Section 3.1 for the 
definition of the excess variance). In such cases 
a meaningful excess variance cannot be obtained 
and its value is consistent with being zero (The 
same situation also happened in so me objects or 
observations in previou s studies, e.g. lO’Neill et all 
2005 ; Ponti etahli2012l) . We find that such sources 


2 Except the object J1153+4612, since it is bright enough 
(with a mean count rate=2.1 counts s —1 ) that the afore¬ 
mentioned influence on timing analysis can be ignored. 
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Fig. 1.— Example light curves of low-mass AGN 
in our sample observed with XMM-newton (panels 
1-4) and ROSAT (panel 5), respectively. All the 
light curves are extracted with 250 s time bins. 


mostly have the mean S/N0 below 2.3. We thus 
introduce a cutoff on the mean S/N, below which 
the sources are dropped 0 . Our final sample af¬ 
ter the S/N cut is composed of 11 objects with 15 
observations. Among them, 10 objects were ob¬ 
served with XMM-Newton for a total of 13 obser¬ 
vations, and 2 objects observed with ROSAT for 
2 observations (J1223+0726 was observed both in 
XMM-Newton and ROSAT observations). Table 
□ summarizes the basic parameters of the sample 
sources and the information on the X-ray obser¬ 
vations. The black hole masses are taken from 
, which are in the 10 5 — 10 6 Mq 
range and the accretion rates in the 0.06 — 0.90 
range (Figure [2]). All the objects are all at very 
low redshifts 2 < 0.21 with a median z = 0.090. 


Dong et ajj (2012) 



log(M B H/M G ) logtLaa/LEj,) 


Fig. 2.- Distributes of the black hole mass and 
sca led accre tion rate for our sample, taken from 
Dong et al. ( 2012 1. 


3. MEASUREMENT OF EXCESS VARI¬ 
ANCE AND THE PSD MODELS 


3.1. Excess Variance and the Uncertainty 


FollowinglNandra et al. (199711 (see alsolTurner et al 


19991: Vaughan et al.ll2003a: Ponti et al.ll2004ll . the 

normalized excess variance is calculated using the 
definition, 




1 N 
AT/? II 


[(X z - fi) 2 - af 


(1) 


3 The mean S/N is calculated as the ratio of the mean source 
count rate to its mean statistical error. 

4 One object with a signal to noise ratio (S/N~2.5) above the 
threshold (J1720+5748, observed with ROSAT) turns out 
to have the intrinsic variability consistent with zero. This 
is likely due to the stochastic nature of the X-ray variability 
of AGN, which happens to result in little variations within 
a relatively short time span. This object is also excluded. 
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Table 1 

Sample and information on X-ray observations 


Num 

(1) 

SDSS Name 

(2) 

z 

(3) 

log(M B H) 

(M 0 ) 

(4) 

Lhol/ -^Edd 

(5) 

ObsID/SEQID 

(6) 

Count rate 
(counts s — 1 ) 

(7) 

Expo. 

(ks) 

(8) 

2 

°rms,10ks 

(9) 

FYlTOks 

(10) 

XI 

J102348.44+040553.7 

0.099 

5.44 

0.32 

0108670101 

0.08 

51.00 

0.044 

0.019/0.019 






0605540201 

0.13 

106.50 

0.057 

0.012/0.012 

X2 

.1114008.72+030711.4 

0.081 

5.70 

0.90 

0305920201 

0.75 

39.00 

0.015 

0.004/0.004 

X3 

J122349.64+072657.9 

0.075 

5.63 

0.55 

0205090101 

0.19 

24.00 

0.027 

0.025/0.025 

X4 

J143450.63+033842.6 

0.028 

5.73 

0.06 

0305920401 

0.12 

22.00 

0.101 

0.048/0.048 






0674810501 

0.13 

11.75 

0.065 

0.032/0.032 

X5 

J010712.03+140845.0 

0.077 

6.09 

0.34 

0305920101 

0.31 

19.00 

0.055 

0.027/0.027 

X6 

J135724.53+652505.9 

0.106 

6.20 

0.47 

0305920301 

0.51 

18.75 

0.040 

0.017/0.017 






0305920601 

0.49 

13.75 

0.014 

0.007/0.007 

X7 

J082433.33+380013.2 

0.103 

6.11 

0.41 

0403760201 

0.09 

15.00 

0.097 

0.059/0.058 

X8 

J091449.06+085321.1 

0.140 

6.28 

0.37 

0744220701 

0.76 

31.75 

0.040 

0.011/0.011 

X9 

J134738.24+474301.9 

0.064 

5.63 

0.67 

0744220801 

0.63 

20.75 

0.028 

0.009/0.009 

X10 

J115341.78+461242.3 

0.025 

6.09 

0.29 

0744220301 

2.16 

13.75 

0.019 

0.007/0.007 

R1 

J122349.64+072657.9 

0.075 

5.63 

0.55 

RP600009N00 

0.10 

16.00 

0.186 

0.077/0.076 

R2 

J111644.65+402635.6 

0.202 

5.83 

0.45 

RP700855N00 

0.08 

18.50 

0.032 

0.030/0.029 


Note. —Column 1: X or R denotes the object observed with XMM-newton or ROSAT , respectively, and X9, X10, XI1 are new observations 
from our programme (proposal: 074422 , PI: W. Yuan); Column 2: object name; Column 3: redshift; Column 4: black hole mass in units of the 
solar mass Mq, from [Dong et ahl d2012f) ; Column 5: Eddington ratios; Column 6: observation ID for XMM-Newton or sequence ID for ROSAT 
observation; Column 7: mean count rate of each observation (countss —1 ); Column 8: cleaned exposure time (ks); Column 9 &; 10: excess variances 
and errors calculated on a timescale of 10 ks. 









where N is the number of good time bins of an X- 
ray light curve, fi the unweighted arithmetic mean 
of the count rates, X t and er* the count rates and 
their uncertainties, respectively, in each bin. 


As shown bv lvan der Klisl ( 19891 1997 ). the ex¬ 
cess variance is the integral of the PSD of a light 
curve over a frequency interval given by Eq. (2), 


cr 


2 

rms 


P (is) dis, 


( 2 ) 


where is m i n = y, v max = 2 AT’ ^ an( 4 ^T are the 
time length and the binsize of the light curve, re¬ 
spectively. For a given light curve, it is clear that 
the exact value of the excess variance is depen¬ 
dent on the length of the light curve (e.g. obser¬ 
vational duration) as well as on the binsize AT. 
In order to compare the excess variances of dif¬ 
ferent objects or observations in our sample, the 
duration (timescale) and the binsize of the light 
curves should be set to be the same for all the 
obejects. For this purpose, the light curves are di¬ 
vided into one or more segments of 10 ks in length, 
and for each of the segment the excess variance is 
calculated. For observations having more than one 
segment, the mean of the segments is taken. The 
results are listed in Tabic [l] 


The uncertainty of the excess variance comes 
from two sources, one of the measurement un¬ 
certainty and the other of the stochastic na- 
ture of the v a riabilit y process, as shown by 
Vaughan et ali (;2003ah . The measurement un¬ 


certaintie s of the excess varian ce are estimated 
following 1 Vaughan et al. (2003a): 



(Ac 


where (erf') is the mean of the square of count 
rate uncertainties, F var the fractional variability 
(F var = of ms ), and the other quantities are de¬ 
fined in Eq. (1). 

As discussed in Vaughan et al. (l2003al) . the 
light curves of AGN are simply stochastic, mean¬ 
ing that each observed light curve is only one real¬ 
ization of the underlying random variability pro¬ 
cess, and each realization can exhibit a slightly dif¬ 
ferent mean count rate and variance. The random 
fluctuations between different realizations lead to 
non-ignorable scatter in the excess variance. This 


phenomenon can be reflected in our result: for 
a source with more than one observation (such 
as J1023+0405, No. XI in Table [U, the excess 
variance varies considerably. The scatter can be 
significantly reduced if the light curves are suffi¬ 
ciently long, or, equivalently, having a sufficient 
number of data points. We then compute the un¬ 
certainty owing to the stochastic nature of the 
va riability proces s using the method introduced 


Vaughan et al. (12003aD . We first build a PSD 


m odelyl to simulate a light c urve using the method 
of lTimmer fc Koenig! (1 9951) with a binsize of 250 s 
and sufficiently long duration. The light curve is 
divided into 1000 separate segments of 10 ks. The 
distribution of the calculated excess variances of 
all the segments is obtained, from which the range 
of stochastic uncertainty is found. We take the 
68% confidence range for the expected excess vari¬ 
ance on each of the timescales to get the stochas¬ 
tic scatter. Thus, from our simulations the scatter 
due to this random process is A log(of ms ) = —0.29 
and +0.24 for 10 ks light curve of a 250s-binsize. 


Finally, the total uncertainties of the excess 
variance are obtained by combining in quadra¬ 
ture the stochastic and measurement uncertain¬ 
ties, which are at the 68% confidence level. The 
obtained uncertainties of of ms are given in TablefT] 
For an object having more than one observation, 
the mean and its uncertainty are calculated and 
used as the final excess variance in the following 
analysis. 


3.2. PSD Models of AGN 

As mentioned above, the excess variance is the 
integral of the PSD of a light curve over a fre¬ 
quency interval. If the shape of the PSD is known, 
it is possible to calculate the values of the excess 
variances, and thus to derive the model relation 
between the excess variance and black hole mass, 
which can be compared to observations. For AGN, 
it has been shown that the PSD function has a 
standard shape of a broken power-law, 


P(u) 


A(v/v br) \ V < Z+r 
A(v/v br)“ 2 , V > Z+r 


(4) 


5 The PSD models used in this paper will be introduced in 
the next subsection, and all the models yield essentially the 
same results. 
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where ty is the break frequency (iMcHardv et al 


2006; Gonzalez-Martm fo Vaughan 20 1 2t) . For an 
AGN with given Mbh and to, zy can be deter¬ 
mined, though the exact value differs somewhat 
in different models. From Eq. (2) and (4) the ex¬ 
cess variance can be expressed as 


f'l ^br(A... Vyn -L 


c 1 

Clin 


), 


In (-at*-) - ( -22a-') + 

y ^min J Y I'max J 

( Cm ax 1 

Cm in J 5 


^br A Z-bniri 
Z-'min < Z^br < 
Z-V > ^max 

(5) 


servational data, and a new to dependent 
PSD amplitude was suggested: C'i = am ~P, 
where a = 0.003±g;ggi and /3 = 0.8 ± 0.15 
were fitted in their study. 

Model C: Similar to model A (Ci = 0.017), 
but an improved scaling relation v b r = 
0.001 to°- 24 (Mrh / 10 6 AC ,) -1 suggested by 


Gonzalez-Martm & Vaughan (2012), is adopted 


where Ci = Azy is defined as the PSD am- 


plit u de (Panadakid 20041 s ee also model A in 
Gonzalez-Martm et al. 20 111 for details). It has 
been suggested that tv is inversely scaled with 
Mbh, with a possible dependence on the scaled 
accretion rate in the Eddington units. The rela¬ 
tionship between PSD and of ms is illustrated in the 
sketch in Figure [3J in which the PSD of three ob¬ 
jects in our sample or P12 are shown. The break 
frequencies derived from Mb h and to using the 
Gonzalez-Martm & Vaughan ( 20121) scaling rela¬ 
tion are also indicated for the three objects, which 
represent the three cases in Eq. (5), respectively. 

However, although the broken power-law model 
given in Eq. (4) is widely accepted, the exact 
values of the amplitude C\ and the break fre¬ 
quency zy differ som ewhat in different stud- 


Papadakisl ( 200 4) first suggested a model 


les. 

with zy = 17/(Mbh/M 0 ) (Hz) and Ci = 

Av br = 0.017. Later studies with larger sam¬ 
ples found a secondary de pendence of zy on the 
scale d accretion ra te m (IMcHardv et al.l 120061 : 
Gonzalez-Martm fc Vaughan! 12012 ). in addition 
to the dependence on Mbh- In the literatures 
of recent studies, three following PSD models of 
AGN were used: 


Model A: Adopted by iGonzalez-Martfn et al.l 
( 201ll) . this model assumes the scaling 
relation zy = 0.003 to(Mrh / 10 6 M ^) -1 
suggested by IMcHardv et all (i2006i ). and 
a constant PSD amplit ude C'i = Azy 
(= 0.017 ± 0.006. |Papadakisll2004lh 


Model B: The scaling relation for zy is the same 


as in model A. As argued by iPonti et al 


( 20121) . the dispersion of the Mbh — of ms re¬ 


lation arising from different to derived from 
model A is too large to account for the ob¬ 


ax which predicts a weaker dependence on to. 

4. RESULTS AND DISCUSSION 

4.1. Mbh — of ms Relation in the Low-mass 
Regime 

The measured excess variances of 11 low-mass 
AGN of our sample are calculated from the X-ray 
light curves obtained in Section 3 using Eq. (1). 
The results are listed in Table |T| To enlarge the 
sample, we make use of both the X'.MM-Newton 
and ROSAT data, which have somewhat differ¬ 
ent energy bands. It has been shown that the 



calculated using a binsize of 200 s, different from 
ours (250 s). However, the effect of such a differ¬ 
ence on the calculated excess variance is negligible 
(< 2%), as found by performi ng s imulations using 
the method of lTimmer fc Koenig! ( 1995i ). 


Here we study the relation between cr 2 ms 
and Mbh by combining the low-mass sample 
and the P12 sample. The combined sample 
has Mbh spanning four orders of magnitude, 
10 5 — 10 9 M 0 . The results are shown in Fig¬ 
ure U Also plotted is the relation [log(of ms ) = 
—2.09 — 1.031og(MB H/10 7 Mfb)1 derive d from the 
P12 sample only bv IPonti et al. ( 20121) . It shows 
that the cr 2 ms values of the low-mass sample are 
comparable to the largest values in P12 with 
higher Mbh- However, for all the sources in the 
combined low-mass sample except one, the excess 
variances fall systematically below the extrapola¬ 
tion of the Mbh — cr 2 ms relation derived from the 
high-mass P12 sample. Our result indicates that, 


6 They are GH18, GH49, GH78, GH112, GH116, GH138, 

GH142, and GH211. 
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in the low-mass regime, the Mbh — ff r 2 ms relation 
may deviate from the previously known linear re¬ 
lation, and is likely to flatten at around ~ 1O 6 M 0 
toward the low-mass end. In fact, the low-mass 
sample objects themselves do not show any cor¬ 
relation between of ms and Mbh (the Spearman’s 
correlation test gives a null probability of 0.42). 
To quantify the statistical significance of this de¬ 
viation, we perform two statistical tests, assuming 
that the previous inverse linear relation is a good 
description of the Mbh — q 2 ms relation over the 
entire mass range. First, we use the two-sided 
binomial test to find out the probability of having 
18 out of 19 low-mass AGN fall below the relation 
as observed, whereas an equal probability (50%) 
of falling on either side of the relation is expected. 
This gives a probability of 7 x 10 -5 . Moreover, 
we fit an inverse linear relation with the slope 
fixed at -1 in the log-log space to the data of both 
the low-mass and the P12 sample. The reduced 
X 2 /dof of the fit is 836/53. Then we assume that 
the inverse linear relation breaks to a constant 
at masses below a so-called break mass, which is 
fitted as a free parameter. The y 2 /dof value for 
this model is 705/52 (the break mass is fitted to 
be 1.7 x 10 6 Mq). This improves the fit dramat¬ 
ically by reducing the y 2 by Ay 2 = 131 for one 
additional free parameter. Despite the fact that 
although the fitting is statistically not acceptable 
in terms of the large reduced y 2 (may arise from 
some intrinsic scatters inherent to the relation- 
ship ), as an appro ximation, we employ the F-test 
( Bevineton fc Robinsonl Il992n to test the signif¬ 
icance of adding the flattening term does NOT 
improve the fit. This yields a small p -value 0.003. 
We thus conclude that the flattening of the inverse 
linear Mbh — of ms relation toward low-mass AGN 
is statistically significant. The fitted break mass 
is - 1.7 x 1O 6 M 0 . 


We also fit a simple linear relation with the 
slope as a free parameter in the log-log space 
using the data of both the low-mass and the 


Ponti et al. ( 

2012) 

Bianchi et al. 

2009 


(see also Appendix A of 


over-plotted in Figure [4] as dotted lines. We 
find log(cr r 2 ms ) = (-2.09 ± 0.044) + (-0.58 ± 
0.053) log(M BH /lO 7 M 0 ). The slope -0.58i0.053 
deviates significantly from the expected value -1 
based on previous studies for AGN samples with 



Frequency (Hz) 


Fig. 3.— Sketch of the power specturm density of 
AGN and its relationship with the excess variance, 
which is the integral of the PSD over the frequency 
span [u m in, r'max], where ^ max is set by the binsize 
(250 s) and z/ m ; n is set by the length (10 ks) of the 
light curve. The PSD of three AGN with differ¬ 
ent Mbh and thus differnt break frequency iy, r are 
shown as examples, which correspond to the three 
Cases in Eq. (5) of Z^br ^ Z/ m i n , ^min ^ A ^maxj 
Vbr > Z'max, respectively. 



Fig. 4.— Relationship of the black hole 

mass and the excess variance, which are mea¬ 
sured on a timescale of 10 ks, for our sam- 
ple objects (s t ars), non-overlapping sources from 
Ludlam et al. (12015) (op en circles), and the sam¬ 


ple of Ponti et al. ( 2012 1 (dots). The errors are 
at the 68% confidence level. The solid line is the 
best-fi t relation (slope ss -1) given by Ponti et al 


( 20121 ) based on that sample only, and the dotted 
line represents the best-fit relation (slope ~ -0.58) 
for both the low-mass and the P12 sample (see 
text). 














































highe r Mbh (e.g. Zhou et al. 20ld : Ponti et al 


20121 ) . 


For low-mass AGN, there appears to be a large 
scatter in the of ms values spanning almost one 
decade. The distribution of the excess variances 
(in logarithm) for our sample objects with Mbh < 
2 x 1O 6 M 0 is plotted in Figured] (left panel), along 
with a fitted Gaussian distribution (dashed line). 
It would be interesting to examine the intrinsic 
dispersion of their of ms distribution, by taking 
into account the uncertainty of individual of ms . 
The maximum-likelihoo d method as introduced by 
Maccacaro et al. (1 988f ) is used to quantify the in¬ 


trinsic distribution (assumed to be Gaussian) that 
is disentangled from the uncertainty of each of 
the measured excess variance (also assumed to fol¬ 
low a Gaussian distribution). We find a mean 
(log(af ms )) = —with a standard devi¬ 


ation er = 0.22l o ;gg®. The confidence contours of 
the two parameters are shown in Figure [5] (right 
panel). This indicates that the apparently large 
scatter of of ms can mostly be attributed to the un¬ 
certainty of each of the measurements (including 
both the measurement and the stochastic uncer¬ 
tainties), and the intrinsic dispersion is small, but 
non-negligible (the standard deviation of log(of ms ) 
ss 0.15 — 0.3 at the 68% confidence level). It may 
also suggest that the dependence of of ms on any 
other parameters (e.g. accretion rate) is likely not 
strong. 


In conclusion, the previously found inverse 
scaling of Mbh — of ms for AGN with supermas- 
sive black holes cannot be extrapolated to the 
low-mass regime, but starts to flatten at around 
Mbh ~ 1O 6 M 0 . This implies that, although the 
excess variance can still be used to search for AGN 
with A/bh ^ 10 6 M 0 , it fails to provide accurate 
black hole mass estimation in this mass range. 
This result is co nsistent with t hat ob tained in a 
recent paper by iLudlam et al. (2015!). It is also 
suggested that, for AGN in the 10 5 — 10 6 M 0 
range, the intrinsic dispersion of the excess vari¬ 
ances is likely to be small [standard deviation of 
logf^rms) ~ 0.15 — 0.3]. 


4.2. Comparison with Model Predictions 

Here we compare the observed Mbh — af ms rela¬ 
tion with the predictions based on the above three 
PSD models. In theory, the Mbh — of ms rela¬ 


tion can be deduced from the PSD of AGN using 
Eq. (5), once the shape and parameters of the PSD 
model are known. We take the three PSD models 
in Section 3.2, and fix the Vbr(M bh, to) relation as 
their original forms. Since the PSD amplitude Ci 
was determined by f itting the Mbh — oW,„ relation 
in pr evious work ( Paoadakisl 201)4 Ponti et all 
20121 ). it is set to be a free parameter here (two 
parameters a and /3 for model B since C\ = 
am _/3 ). The bivariate models of ofms(M bh, bi) 
from Eq. (5) are fitted to the data of both the 
low-mass and the P12 samples using the simple 
^-minimization fitting. The results are given in 
Tabic 0 with the errors of the fitted normalization 
are at the 68% confidence level. The best-fit PSD 
amplitudes are very close to their original values 
for all the models. 


The model Mbh — of ms relations using the best- 
fit PSD normalizations are plotted in Figure [6] for 
the three models, along with the data. The black, 
red and blue solid lines represent the relations for 
three typical accretion rates (to = 0.02,0.29,0.97), 
respectively, which are the mean (in logarithm) of 
to in three to bins (0.01-0.1, 0.1-0.5, and 0.5-4.0). 
The objects with m in the three to bins are plotted 
in corresponding colors. It can be seen that all the 
three models can reproduce the observed trend of 
the A/bh — of ms relation well quantitatively. For a 
given accretion rate, an inverse proportion is pre¬ 
dicted in the high A7 bh regime, whereas in the low- 
mass regime it flattens toward lower Mbh- The 
exact value of A7 bh at which the relation starts to 
flatten (the break mass) depends on iin: the higher 
to, the larger the break mass. This trend generally 
holds for all the three PSD models, although the 
break mass varies from model to model, in addi¬ 
tion to its dependence on rh. 

Besides this general trend, there are some no¬ 
ticeable differences between the three models. For 
model A, the strong dependence on to leads to a 
large scatter in high-mass regime over a range of 
to. Model B predicts a smaller scatter, but the 
scatter is relatively larger in the low-mass than 
in the high-mass regime. Compared to model A, 
model C also predicts a smaller scatter over the 
whole Mbh range, given its weak dependence of 
V br on TO. 

The inverse proportion of the Mbh — of ms re¬ 
lation and its flattening toward low A/bh can eas¬ 
ily be understood from Eq. (5) and Figure [3] For 
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Table 2 

Fitting results of the M b h — relation with the three PSD models. 


PSD model 

PSD amplitude 

X 2 /do,f 

(i) 

(2) 

(3) 

A 

Ci = 0.032 ± 0.001 

2048/53 

B 

a = 0.0074 ± 0.0004, /3 = 0.80 ± 0.02 

598/52 

c 

Ci = 0.022 ± 0.001 

820/53 



'ogOD 


0.5 



O 

% 0.3 
"> 

CD 
TJ 

B 0.2 

o 
"a 
c 
o 

o.l 

0.0 

- 1.7 - 1.6 - 1.5 - 1.4 - 1.3 - 1.2 - 1.1 

<log((7r ms )> 


Fig. 5.— Left panel: distributions of the excess variances for the objects with M BB < 2 x 10 6 Mq. The 
dashed line represents a fitted Gaussian with a stardard deviation of 0.28. Right panel: confidence contours 
(at the 68%, 90% and 99% confidence levels) of the mean and standard deviation of the intrinsic distribution 
of (solid lines), which are derived using the maximum-likelihood method (see the text). The cross marks 
the best-estimated values. 




Fig. 6.- The same as Figure [4] except that the solid lines represent the theoretical relations derived from 
the best-fit PSD models assuming different normalization and the dependence of break frequency on black 
hole mass and accretion rate (model A, B, and C; see text for details). The three colors correspond to three 
accretion rates of to =0.02 (black), 0.29 (red), 0.97 (blue), which are chosen to be the mean (in logarithm) 
of 77i in three rh bins 0.01-0.1, 0.1-0.5, and 0.5-4.0 of the whole sample objects. The sources with to in the 
three ifi bins are plotted in corresponding colors. 


10 
































































J'br < v min (e.g. NGC 3227 in Figure[3]), as I'br in¬ 
creases (and Mbh decreases, since i^, r oc Mbh^ 1 ), 
the integral of the PSD within [j/ m in, ^max] in¬ 
creases accordingly, resulting in a linearly increas¬ 
ing cr 2 ms . Thus the previously known Mbh — cr 2 ms 
relation (for Mbh > 10 6 Mq) is a manifestation of 
the more fundamental dependence of vy, v on Mbh 
(see also IPaoadakisl 2004 : Ponti et al. 2012ft . It 
also shows that the secondary dependence of tv 
on to introduces a scatter, the extent of which de¬ 
pends on the models. For v m \ n < v b r < ^ max (e.g. 
MRK 766), the relation becomes non-linear. For 
Vb r > i/max (e.g. J1023+0405), cr 2 * ms becomes in¬ 
dependent of v b r (thus of Mbh) and remains a 
constant. 

We have shown that the observed Mbh — of ms 
relation over almost the whole Mbh range for 
AGN of Mbh = 10 5 — 10 9 M 0 can be explained 
qualitatively by the PSD shape of AGN and the 
dependence of v\,r on Mbh, although the exact re¬ 
lation depends on the details of the models. It 
would be interesting to investigate which of the 
above models give a better description of the data, 
taking advantage of the extended Mbh range pro¬ 
vided by our sample. Here we discuss this issue 
only briefly by comparing the deviations of the 
data from the models in terms of the fitted y 2 , al¬ 
though they are too large for the fits to be accept¬ 
able nominally for all the three models (the large 
y 2 values are partly due to the somewhat large un¬ 
certainties in Mbh and to, which are not taken into 
account in the fitting). We discuss the three mod¬ 
els, respectively, (i) Model A has a much larger y 2 
value than model B and C, and is likely not a good 
description of the data. The ame suggestion was 
also argued bv IPonti et al. ( 2012ft . (ii) For model 
B, The fitted param eter 0 (= 0.8 0 ± 0.02) is con¬ 
sistent with that of Ponti et al. ( 2012h . For the 
objects with Mbh < 2 x 10 6 M 0 , no correlation 
is found between cr 2 ms and to (the Spearman cor¬ 
relation test gives a null probability p = 0.42), 
which seems to be inconsistent to the model pre¬ 
diction (see Figure [H middle panel). However, 
this may be partly due to the uncertainties in 
determining rh. (iii) For model C, our result 
yields a P SD amplitude consistent with that of 
IPapadakfsl ( 2004h . which is only weakly dependent 
on to. This leads the excess variances converg¬ 
ing to a constant value toward low Mbh, which is 
roughly consistent with the small intrinsic scatter 


(not zero, however) found above. Overall, based 
on the fitted y 2 values (Table [2|, we tend to con¬ 
sider models B, and perhaps model C to be better 
descriptions of the PSD shape of AGN. However, 
a rigorous comparison of the PSD models is be¬ 
yond the scope of this paper and will be carried 
out in future work with a larger sample and better 
quality of data. 


5. CONCLUSION 


The Mbh — cr 2 ms relation for AGN in the low- 
mass regime (Mbh < 10 6 M 0 ) is investigated us¬ 
ing a sample of 11 low-mass AGN observed by 
XMM-Newton and ROSAT , including both new 
obser vations and ar c hival d ata, as well as 8 sources 
from Ludlam et al. ( 2015lh We find that, the in¬ 
verse linear Mbh — cr 2 ms relation established in 
the high-mass regime in previou s studies (e.g. 


Zhou et al. 201ft Ponti et al. 2012 1 fails to extend 
to Mbh below 10 6 M 0 . The relation becomes to 
flatten at ~ 10 6 M 0 , below which the excess vari¬ 
ances seem to remain constant. Our result is in 
good agreement with that obta i ned from a recent 
similar study by iLudlam et al. ( 2015 ). This is in 
fact consistent with the model prediction from our 
current understanding of the PSD of AGN and the 
dependence of the break frequency on Mbh- 
Our result suggests that while the X-ray excess 
variance may still be used to search for low-mass 
AGN candidates 0, it fails to provide reliable es¬ 
timation of the black hole mass for AGN with 
Mbh 1$ 10 6 M 0 . In this Mbh regime, it is also 
found that the excess variances show small intrin¬ 
sic dispersion, when their uncertainties (both mea¬ 
surement and stochastic) are taken into account. 


This work is supported by the National Natural 
Science Foundation of China (Grant No.11473035, 
Grant No. 11033007 and Grant No. 11173029), 
and the Strategic Priority Research Program 
“The Emergence of Cosmological Structures” of 
the Chinese Academy of Sciences (Grant No. 
XDB09000000). This work is mainly based on 
observations obtained with XMM-Newton, an 
ESA science mission with instruments and con¬ 
tributions directly funded by ESA Member States 

7 As t hose already carri ed out by e.g. iKamizasa et al ]|m2) 

and (Terashima ( 2 0l3) by assuming that the reverse linear 

scaling relation can be extended to the low-mass regime 
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